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.2 ' Abstract 

Rotation of a single bacterial flagellar motor is powered by multiple stators tethered to the 
cell wall. In a "power-stroke" model the observed independence of the speed at low load on 

> 

' the number of stators is explained by a torque-dependent stepping mechanism independent 

m ' 

of the strength of the stator tethering spring. On the other hand, in models that depend 
qq , solely on the stator spring to explain the observed behavior, exceedingly small stator spring 
constants are required. To study the dynamics of the motor driven by external forces (such 
as those exerted by an optical tweezer), back-stepping is introduced when stators are driven 
far out of equilibrium. Our model with back-stepping reproduces the observed absence of a 
barrier to backward rotation, as well the behaviors in the high-speed negative-torque regime. 
Recently measured temperature dependence of the motor speed near zero load (Yuan & Berg 
2010 Biophys J) is explained quantitatively by the thermally activated stepping rates in our 
model. Finally, we suggest that the general mechanical properties of all molecular motors 
(linear and rotary), characterized by their force(torque)-speed curve, can be determined by 
their power-stroke potentials and the dependence of the stepping rates on the mechanical 
state of the motor (force or speed). The torque-speed curve for the clockwise rotating 
flagellar motor has been observed for the first time recently (Yuan et al. 2010 PNAS). Its 
quasi-linear behavior is quantitatively reproduced by our model. In particular, we show 
that concave and convex shapes of the torque-speed curve can be achieved by changing 
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the interaction potential from linear to quadratic form. We also show that reversing the 
stepping rate dependence on force (torque) can lead to non-monotonicity in the speed-load 
dependency. 



INTRODUCTION 

Like all other motile bacteria Escherichia coli continuously explores its environment in 
pursuit of favorable conditions for its survival. E. coli cell's swimming motion is propelled 
by the rotations of several flagella, each driven by a rotary flagellar motor attached to the 
cell body. The cell changes its direction of motion by switching the rotational direction of 



its flagellar motors jl-|3j. When all flagellar filaments are driven to rotate counterclockwise 
(CCW), they form a bundle that propels the cell forward in a smooth swimming ("run") 
mode. When at least one of the flagella rotates clockwise (CW), the coherent bundle disas- 
sembles and the cell tumbles as a result. The cell's signaling pathway enables it to measure 
and compare environmental changes over time. Positive changes -increases of chemoattrac- 
tant - prolong the duration of the smooth swimming motions, so that the random walk of 
the cell is biased towards a preferred direction, i.e., the direction with higher attractant 
concentrations. 

Each flagellar filament (~ lO/xm long) is driven by a rotatory motor embedded in the 
cell wall, with an angular speed of the order of 100 Hz A single motor is composed of 
one rotor, a circular ring structure roughly 45nm in diameter 4|, and several stator units 
anchored to the rigid peptidoglycan cell wall. The rotation of the rotor drives the flagel- 
lar filament through a short hook. The rotor contains a ring of ~ 26 FliG proteins and 
each stator has four copies of proteins Mot A and two copies of proteins MotB, forming 
two proton-conducting transmembrane channels. A flow of protons, driven by the electro- 
chemical gradients across the channels, causes conformational changes of the stator proteins 
that generate force on the rotor through electrostatic interaction between MotA and pro- 
tein FliG j^]. In some alkalophilic and marine Vibrio species of bacteria, sodium ions take 
the place of protons. The work per unit charge that a proton can do in crossing the cyto- 
plasmic membrane through the proton channel is called the "proton- motive force" (pmf). 
Which consists of two terms: a transmembrane voltage Aip and an entropic contribution 
due to the ion concentration difference across the cytoplasmic membrane proportional to 
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ApH=log( [H+ riplasm ] / [H+ ytoplasm ] ) : 

Ail = PMF = Aip + 2.303— ApH. (1) 

e 

In the case of a sodium driven motor, ApH is replaced by the sodium ion concentration 
term ApNa = -/o5f([iya + ]cytopiasm/[A^a + ]p ei . iplasm ). 

The duty ratio of the flagellar motor is close to unity [61], which means that at any given 
time, at least one stator is engaged with the rotor. Biochemical and structural studies indi- 
cate that the motor torque is generated by stator conformational changes as a consequence 
of ion binding and unbinding to the negatively charged D32 residue on the MotB helices 
(D24 on PomB for the Na + motor). Only few portions of the atomic structure of the rotor 
are known and details of the stator-rotor interaction remain unclear due to lack of 

structural information. Presumably, the passage of (one or a few) protons switches the stator 
to be engaged with the next FliG monomer on the FliG ring along the direction of rotation, 
stretching the link between the stator and the rotor. The subsequent relaxation process 
rotates the rotor and the attached load towards the new equilibrium position. This can give 
rise to a step-like motion, characterized by advances of the rotor followed by waiting periods, 
as has been demonstrated by direct observation [lOj for a sodium-powered motor at very 
low "sodium-motive-force" . However, a general understanding of the stepping dynamics of 
a single flagellar motor is still lacking. 

The torque-speed dependence is the key biophysical property that characterizes the motor 



12jj. For E. coli at room temperature under physiologically relevant conditions, the 
maximum angular velocity is ~ 300Hz and the estimated maximum torque ranges from 

nn 

1400 pN-nm [13j, [14| to 4700 pN-nm [151 ] . The measured torque-speed curve for bacterial 
flagellar motor in the CCW rotational state shows two distinctive regimes: i) a plateau 
at high load-low speed and ii) a fast quasi-linear decay at low load-high speed. At the 
transition between the two regimes the angular velocity reaches a large fraction (« 60%) of 
the maximum velocity, the torque falls by only ~ 10% of its maximum value, r max , and the 
torque-speed curve exhibits a characteristic concave shape, like a "knee" with the leg pulling 
forward. The general mechanism responsible for the shape (concavity) of the torque-speed 
curve remains to be explored, particularly in connection with other molecular motors, such 
as kinesin and myosin-V, where a convex, at high load, and a concave, at low load, force- 
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for the CW 



speed curve was observed [lg, |l7|. In addition, new experimental results 
rotation case show a different torque-speed curve shape respect to the CCW case. 

Another important characteristic of the torque-speed curve is the absence of a barrier to 



backward rotations. Electrorotation 



191 ] and optical tweezer [15| experiments show almost 



a constant value for the torque near zero speed, when the motor is driven backwards. The 
torque does not increase and remains at a constant value in a range up to 100Hz in both 
directions of the rotational speed. The mechanism for the absence of a backward rotation 
barrier is unclear. 

Several mathematical models 20M25| have been proposed to explain various aspects of the 



observed torque-speed characteristics based on assumptions about details of the electrostatic 

42 9| use a general 



interaction between the stators and the rotor. More recent works [26 
approach independent of the details of the microscopic interactions. The model introduced 
in [26| can reproduce the observed torque-speed curve characteristics. However, it predicts 



the maximum velocities depend inversely on the number of stators, Supporting Text in 
which is inconsistent with the observed behavior at low load, where a recent experiment 
shows that near zero load the velocity of the motor is independent of the number of stators. 
Recently, we have introduced a model of the motor in which the stepping probability of 
a stator depends on its interaction with the rotor [31 1 - We show that under the general 
assumption that the stepping rate is higher when the stator is generating torque in the 
"wrong" direction, i.e., opposite to the rotational direction of the rotor, the zero-load speed 
at zero load is independent of the number of the stators. An alternative explanation 27, 3o| 
was put forth by introducing a spring between the stator and its anchoring point at the rigid 
cell wall. The idea to model the stator as a spring anchored to the cell wall, was introduced 
for the first time in (32). The flexible stator location enabled by a soft spring can damp the 
conflicting interactions between stators at low load. However, in this explanation is missing 
a precise study of the chemical and mechanical parameter values, such as the spring stiffness, 
in order to verify that not only the zero-load speed independence on the stator number is 
obtained, but also that the motor show a realistic behavior in all of its parts as well that 
these parameters values are reasonable and cover a sufficient large interval. 

Here, we study the effects of the stator spring systematically. We show that the presence 
of the stator spring does not interfere with the torque-dependent stepping mechanism. On 
the other hand, in models that depend solely on the the soft stator spring to explain the 
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observed behavior, only a narrow range of small torsional constant spring k$ can reproduce 
the experimental observation with reasonable stator position fluctuation and displacement. 

As described before, the flagellar motor can be driven out of its normal operational 
range by external forces. The measured torque-speed relationship in these forced regime can 
inform us about the nature of the stator-rotor interaction and the stator stepping mechanism 
and test different model predictions. For the negative velocity regime (driven backwards 
by external forces, such as a optical tweezer), the torque-speed curve is observed to be 
continuous without any barrier, which is inconsistent with the predictions from models that 



use periodic potential to describe the stator-rotor interaction [26|, [27|, |29|] (see for instance 
Figure 4B in 26] for the potential shape and Supplementary Information for a detailed 
discussion). Here, we extend our model to introduce the possibility of back-stepping which 
becomes significant as the system is driven far away from the stator-rotor equilibrium. This 
extension allows us to analyze the negative speed regime; and our model reproduces the 
absence of a barrier to backwards rotation. We also explore the negative torque region 
where the motor is driven forward by the external force beyond its intrinsic zero-load speed. 
There, our model again reproduces the experimentally observed behavior 33[, i.e., the torque 
first drops linearly below zero before leveling off at a maximum negative torque, comparable 
in magnitude to the positive maximum torque generated by the stator. 

Besides pmf, the motor dynamics is also affected by the temperature. We use our model to 
study the temperature effects by assuming a thermally activated stepping rate and compare 
the model results with the recently measured temperature dependence of the speed at near 
zero load. The resulting activation energy informs us about the possible gating mechanism 
for the flagellar motor stepping. Finally, our simple model can be used to address the general 
question of how the form of the potential describing the interaction between stator and 
rotor determines the torque-speed curve, in particular its general shape. For example, while 
the quasi-linear rotor-stator potential results to a concave torque-speed curve, a parabolic 
potential gives rise to a torque-speed curves with convex or nearly linear shape, similar to the 



one observed in the CW rotational case 



. Force- velocity curves observed in linear motors 



such as kinesin-1 and myosin- V also show convex shapes at high load and concave shapes 
at low load Q, [ill Q| This general relationship provides us a useful tool to predict, either 
the potential shape when the torque-speed/force- velocity curve is known, or the torque- 
speed/force- velocity curve when the potential shape is known. 
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RESULTS 



The power-stroke model of bacterial flagellar motor 



We use the general model framework in our previous work [3JJ with several important 
new ingredients. We now include the probability of stator back-stepping and the stator 
position is anchored to the immobile cell wall through a spring. We first briefly recount 
the main features of the general model. In Figure 1A, a rotor and a stator are shown 
with the interacting potential. The two springs represent the load-rotor interaction and 
the elastic link between the stator units and the peptidoglycan cell wall. Each stator has 
two force-generating subunits. The interaction potential between the rotor and one of the 



two stator subunits is represented by the red curve. As explained in [3l| the two torque- 



generating subunits of a stator interact with the FliG ring (rotor) in a hand-over-hand 



35 



36|. The 



fashion, analogous to the way kinesin proteins interact with microtubules 
forces between the FliG ring and the stators drive the rotation of the rotor. The forward 
hand-switch transition corresponds to a shift of the potential energy in the direction of the 
motor rotation by an angle 5q and the subsequent motor motion is governed by this new 
potential until the next switch of stator subunit occurs. 

The dynamics of the rotor angle and the load angle L can be described by the following 
Langevin equations: 

= ~h E V W ~ ^ " F (° ~ ^ + V^fUa(t), (2) 

i=l 



= f(9 - e L ) + v^ra/3(t), (3) 

where £r and £i are the drag coefficients for the rotor and the load respectively, and N is the 
total number of stators in the motor. V is the interaction potential between the rotor and a 
stator. V depends on the relative angular coordinates AOi = — Of , where Of is the internal 
coordinate of the stator i, which consists of two parts: Of — Oi + 0^. 0{ is the (physical) 
angular distance of the stator from its stationary anchoring point on the cell wall; 0\ is the 
internal (chemical) coordinate of the stator. The load is coupled to the rotor via a nonlinear 
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spring described by a function F, which can be determined from the hook spring compliance 



measurement of ref. 37]]. The last terms in equations ()2]|3]) are stochastic forces acting on 
the rotor and on the load, with ks the Boltzmann constant, T the absolute temperature. 
a(t) and (3(t) represent two independent white noise fluctuations of unity intensity. 

The dynamics of the physical coordinate of the stator 6^ is governed by a Langevin 
equation: 

dOi d 



is ll = ~W V{d ~ ^ ~ Ksdl + V^ffs^t), (4) 

where £s is the drag coefficients of the stator, assumed to be the same for all stators, k$ 
is torsional spring constant of each stator relative to the rotor axis, and 7i(t) is again an 
independent white noise fluctuation of unity intensity. The internal coordinate of the stator 
is changed by the hand-switching process, each switching event changes 9\ by <5 , which is 
half of the angular periodicity of the FliG ring, 7r/26. A forward step (also called a forward 
jump) increases Q\ by Sq. The probability of the forward step during the time interval At is: 
Pf(9f — > 9% + Sq). The stepping rate Pf is controlled by the external driving force (pmf), but 
it also depends on the mechanical state of the stator characterized by the torque between 
the rotor and the stator, r» = —V'(A9i), which depends on the relative angle A9i\ 

Pf{8 C i ->• n + ^o) = Rf(Ti)At = k f (A9 l )At. (5) 

In general, we assume Rffc) to be a decreasing function of r,, e.g., the forward stepping 
rate is assumed to be higher when the stator generates torque opposite to the rotor rotation 



3l| 



In this paper, we introduce the backward stepping probability P\,{9\ —> 9\ — 8 ), which 
also depends on the relative angle A^: 

P b {9t -> 9\ - 8 ) = R b (r t )At = k b (A9 l )At. (6) 

Backward stepping happens when the stator has overreached on the FliG ring with large 
value of A9i (compared with Sq), P b remains near zero at \A9\i < 5q and increases as A9i 
increases beyond Sq. Although back stepping only occurs relatively rarely under normal 
motor operating conditions as the rotation of the rotor prevents very large values of A9$, it 
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becomes more frequent if the rotor is being driven backwards by external forces, which is 
one of the focuses of our paper. 

For the stator-rotor interaction potential V, we have considered different functional forms, 
starting (for simplicity) with the linear ^/-shaped function: V(A6) = r \A6\, where the 
torque from a single stator is r with its sign depending on whether the stator is pulling 
(A9 < 0) or dragging {AO > 0). Correspondingly, the stator forward stepping rate depends 
on the sign of the force: k f (A6 < -5 C ) = 0, k f (-5 c < AO < 0) = k+, k f (A6 > 0) = £;_(> 
k + ), with a cutoff angle 6 C introduced to prevent run-away stators. The form of the backward 
stepping rate will be given in the following sections. Quantitatively, we use Tq = 505pN-nm, 
fa = 0.02pN-nm-s-rad- 2 , fa = 0.004pN-nm-s-rad- 2 , k + = 12000s" 1 , k_ = 2k + , 8 C = S in 
this paper unless otherwise stated. The load fa varies from 0.02 — 50pN-nm-s-rad _1 . Other 
forms of V, such as semi-parabolic and pure parabolic potentials will also be studied, with 
the emphasis on understanding the general relationship between the form of V and the 
corresponding torque-speed curve. 



Independence of the maximum speed on the stator number and the effects of stator 
spring 



In our previous work 31(, it was shown that the maximum rotation speed (at near zero 
load) is independent of the number of stators under the general condition that the stators are 
more likely to step forward when they generate backward torque, i.e., r(= k+jkJ) < 1. We 
call this scenario the force-dependent-stepping (FDS) mechanism. The (physical) locations 
of the stators are fixed in our previous model for simplicity. A more realistic description 
should take into account that the stator is linked to the rigid peptidoglycan through a a-helix 



30. 



32 1 . Since the elasticity of a- helix is known and 



of 7-8nm in length as first suggested in 
its persistence length is about lOOnm [ 



the resulting stator spring constant can be then 
estimated to be at least 500pN-nm-rad _1 , depending on the orientation of the a-helix with 
respect to the rotor. In 27|, it was shown that the independence of the zero-load speed on 
the stator number can be obtained by introducing a stator spring with an elastic constant 
value of k,s ~ 200pN-nm-rad~ 1 or smaller, which is more than 2-fold softer than a 7-8nm 
a-helix spring. In Figure IB it is shown the dynamics of the rotor and the stator under 
large load (Supplementary Figure 1 shows also the case at low load) with such a small stator 
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spring constant as used in [27]. The angular displacement of the stator from its anchoring 
point reach as large as 105, which is 50nm assuming a radius for the rotor of 20nm; and the 
angular position fluctuation reaches 25. The large average stator position displacement and 
its large fluctuation put in doubt this weak-stator-spring (WSS) scenario. 

In this section, we systematically investigate the motor dynamics for a wide range of stator 
spring constants. We emphasize on studying the interplay between the force-dependent- 
stepping and the weak-stator-spring scenarios. We show how adding the stator spring does 
not interfere with the FDS mechanism (r < 1) for a wide range of constant spring values 
considered. However, very stringent constrains for the stator spring value are required in 
the absence of the FDS mechanism, i.e., when r > 1. Such stringent constrains on the stator 
spring constant result from two opposing requirements that need to be satisfied by the 
system. The first requirement is that the stator displacement and the related fluctuation 
should realistically be smaller than a few 5 jsoj, which sets a lower limit for the spring 
constant, k™ m , independent on r. The second requirement is the independence of the zero- 
load speed on the stator number, which requires the spring to be soft and therefore sets a 
maximum value for the spring constant, k™ ax , which depends on r. For r > 1, these two 
requirements select a vanishingly small region of possible k$ values. 

To characterize the dependence of the zero-load speed on the number of stators, we define 
the zero-load speed ratio u r = u; (8)/o;o(l), and the zero-load speed variation e r = |1 — u r \, 
where ojq(1) and o>o(8) are the speed at near-zero load for number of stator N — 1 and N = 8 
respectively. In this case uo(l) = (l±e r )-wo(8). From the recent resurrection experiments at 
low load [30j, the overall variation of the measured maximum speeds for all stator numbers 
is e r ~ 5% in single cell, see for example Figure IB in 30] . Considering the entire range of 
speeds variation in a cell population in the same experiment, a conservative estimate for e r is 
roughly 10%, Figure 2A in [301 ] . In the presence of the FDS mechanism, the zero-load speed 
variation e r is less than 20% for all values of the stator spring constant spring considered 
when r < 1 (Figure 1C), and the variation e r quickly decreases as FDS effects become 
stronger, and when r < 0.5, e r becomes smaller than 10%. If instead it is less probable 
for a stator to step (switch-hand) when it is generating the negative torque, i.e. r > 1, to 
which the model studied in [27] belongs, the zero-load speed variation depends on the stator 
spring constant ks and only becomes small for small Ks- Consider the case r=2, k$ should 
be smaller than 20 pN-nm-rad -1 to have independence of the zero-load speed on the stator 
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number within 10%, and equal or smaller than 200 pN-nm-rad -1 to have independence 
within 15% (three times the value estimated from single cell experiment). In addition, 
soft stator spring leads to large stator position displacement and fluctuations. Figure ID 
shows the stator displacement A8 S between a stator and its anchor point, and fluctuations 
interval Var(# 5 ) around the equilibrium position of the stator (during the waiting phase) as 
a function of the spring constant at high load for the case N=l and r=2 at high load. Values 
less of two 5 correspond to values of Ks bigger than SOOpN-nm-rad^ 1 . On the other side 
bj T inside a zero-load speed variation of 10% needs k,s to be smaller than 20pN-nm-rad _1 , 
which makes a strong argument against WSS's models. All simulation are done with a value 
of the stator diffusion constant of D$ = 500ra<i 2 s _1 . Other values considered of 50rad 2 s~ l 
and 5000ra<i 2 s _1 show very similar behavior. For instance, Supplementary Figures S1A and 
B show the ratio u r as a function of ks for different values of r and with stator diffusion 
constant Ds = 5000ra<i 2 s _1 and D$ = bQradPs^ 1 . 

A more systematic study of the dependence of the model behavior on the two key pa- 
rameters r and k$ is summarized in Figure 2. The FDS condition r < 1 guarantees the 
independence of the zero-load speed on the stator number (within the experimental error) 
for a wide range of stator spring constants. If we add in the criterion that the stator fluc- 
tuation cannot be bigger than the available space for stator, defined as half of the average 
spacing between neighboring stators, which is 22.5° for N = 8, we obtain a conservative 



estimation 



5l| for the minimum spring constant of ~ 488pN-nm-rad l . This value is similar 



to the theoretical estimation of the minimum spring constant for the a-helix. 

On the other hand, in the absence of the FDS condition, only a small region of the 
parameters (r, Ks) that can lead to both small e r and reasonable values of stator displacement 
and fluctuation. The region of acceptable k,s values shrinks as r increases and disappears 
after r becomes bigger than ~ 1.5 (see Figure 2). As well, increasing k$, the initial interval 
of acceptable r between 1 and 1.5 decreases fast and disappears before k s reaches the value 
of ~ S^OOpN-nm-rad- 1 . 



Rotation under external force and the stator's backward steps 

An important characteristic of the torque-speed curve for the flagellar motor is the absence 
of a barrier to backward rotation. Externally applied torques to drive the motor backwards 
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(at negative speed), or forwards at speeds greater than the zero- load speed, have been used 
since the beginnings of the 90's to discriminate between different models. In particular 



thermal ratchet models such the one proposed by Meister 23| predicted that the torque 



will rise sharply if the motor is driven backward. This behavior was first observed [39] but 



successively it turned out to be an artefact of the experimental procedure 40j. Meanwhile, 



Washizu at al. 
either direction. 



4l| found that motor torque is constant up to about 100Hz for rotation in 



In our previous study [3l| we did not consider backward stator steps, which are relatively 



rare under normal operating conditions 10J. However, we expect the back-jumps to become 
dominant when the stator is driven backwards into regimes with A6 < — S c , where the 
forward jumps are prohibited. Here, as described in the Model section of the paper, we study 
the motor behavior under external driving force by including the backward steps (jumps). 
Since the landing points of these back-jumps are still on the positive side of the potential 
with positive torque ro, inclusion of back-jumps in our model can naturally explain the 
actually observed torque continuity near stall when the motor is driven backwards. After 
the first electro-rotation experiment of Washizu et al. [41[ this behavior was confirmed 
using optical tweezer, where the torque generated by the motor was the same when made to 



rotate slowly backward as when allowed to rotate slowly forwards at speeds up to 0.3Hz 15]. 



191 ] showed that there is not barrier to 



Successive and further electro-rotation experiments 
backward rotation at a speed up to 40Hz. In one case, data from one particular cell indicate 
that the linear range for the torque-speed curve extends up to > 100Hz in either directions, 



confirming again the result found in ref. 41]. In both experiments cells were tethered to 
glass coverslip by single flagellum. The E. coli strain used carries a che Y deletion and thus 
rotates flagella exclusively CCW. Such a characteristic of the flagellar motor torque-speed 
curve has been predicted in simple kinetic analysis, involving only stepping rates modeling 
[l^ |. but the study of the negative speed regime has not been addressed so far in more 
general descriptions. Here we explain and interpret these results in the framework of our 
simple model. 

The equation of motion for the load is now modified by an additional term representing 
the external torque applied directly to the load, for instance by mean of an optical tweezer, 



dO 



L 



dt 



F(6 - B L ) - r ext {t) + y/2k B T£ L p(t), (7) 
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where T ext is the constant external applied torque. When a cell is tethered to a surface by a 
single flagellar filament, the motor turns the cell body alternately CW and CCW at a speed 
around 10Hz 42] . This is the starting working point of our simulations, i.e. the value of 
the drag coefficient £x corresponding to this value of the speed, roughly 5pN-nm-s-rad~ 2 , 
is kept constant. Then the torque-speed curve is obtained by steadily increasing r ext by a 
fixed amount for each point of the curve, mimicking the increasing of the electric field in 
electro-rotation experiments as described in At a given value of r ext the value of the 
motor torque is calculated by applying the following torque balance equation: 

£ L UJ = T ext + T. (8) 



For the jumping rate we use the following expression for the forward jump of each stator i: 
k f (A6i) = k g + k af [exp{{ n 56i - E )/k B T)/(l + exp((r i «J0, - E )/k B T))] (9) 

and for the backward steps (jumps): 

h(A0i) = ko^expi-inM - E )/k B T)/(l + ex P {-{r H 5Bi - E )/k B T))], (10) 



where Tj = —V'(A6i) for each stator i, and the force represents the slope of a potential 
equal to the stator- rotor interaction potential but shifted by 0.755 respect to the stator 
position in direction of negative A9i values, see the inset in Figure 3B. 56i and Eq set the 
characteristic torque r c = E /58i, k b = k r kof, with k r a numerical pre-factor, and k g a base 
constant forward rate. As semi-parabolic potential is used here, see Figure 3B. The result 
does not change when a V-shaped potential is used. 

Figure 3B shows the torque-speed curve with N=l when an external torque is applied to 
make the rotor rotate backward and forward. In agreement with the experiments, the motor 
torque varies linearly with speed up to 100Hz in both direction without any barrier (See 
Supplementary Information for the case N=8). It is worth mentioning that including the 
backward jumping probability does not change the qualitative shape and characteristic of 
the torque-speed curve in its normal operation range in which back stepping rarely occurs. 
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In particular, the independence of the zero-load speed on the number of stators is still hold 
(see Supplementary Figure S5 for details). 

Another interesting region outside of the normal motor operating range is the super high 
speed negative torque region when the external torque is applied along the same direction as 
the motor's natural rotation. In this super high speed regime, as the motor speed increases 
beyond its native maximum speed, the torque generated by the motor drops below zero 
linearly with the speed before it levels off at a maximum negative torque as measured by 
|. Figure 3D shows that our model reproduces this behavior. The reason 



Turner et al 



why the torque levels off at a negative torque value equal in absolute value to r can be 
understood from our model as shown in Figure 3C. As the external torque forces the motor 
move faster than its natural maximum speed at low load, the stators spend most of their 
time in the negative torque regime of the rotor-stator potential, where the maximum value 
of the torque is equal to — r . To the extent that the observed maximum negative torque at 
super high speed is comparable to the motor's natural maximum positive torque, our model 
shows that the rotor-stator potential is approximately symmetric. 

The temperature effect at low loads and the thermally activated stepping 

Changing temperature affects all the chemical transition rates, in particular the stepping 
rates. The temperature dependence of k + and in our model leads to changes in the knee 
speed u n and the speed near zero load u without changing the maximum torque at stall, 



which are consistent with previous experimental observations |39l . |43j . In our previous work 
we did not investigate the thermal effects in details. Here, we carry out the study of the 
temperature effects on the torque-speed curve. The main motivation for such a study is given 



by a recent experiments 44] , which extended the previous studies 32j of thermal effects on 
motor rotation at the high-load regime to the near zero load regime. They used nano-gold 
spheres attached to the hooks of the flagella in mutant E. coli cells lacking flagellar filaments 
to measure speed variations of the motor near zero load within a temperature range from 
9 to 37°C. They found that the speed changes nearly exponentially with the temperature, 
and a value of 52 kJ/mol for the activation enthalpy was extracted from the temperature 
dependence of the speeds near zero load. 

For simplicity we consider a V-shaped stator-rotor interaction potential, with the tem- 
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perature dependence of the stepping rates following the simple Arrhenius law: 



k± = k ± exp(-AG/RT), 



(11) 



where R=8.3 JK~ 1 mol~ 1 is the universal gas constant, AG = 52 kJ/mol is the activation 
enthalpy, ko- = 1.6 x 10 13 s _1 , and ko + /ko- = 0.5. Figure 4A shows the dependence of speed 
near zero load on temperature obtained from our model. The agreement with the experimen- 
tal data is excellent with AG given by the experiments and only one fitting parameter ko-- 
At high load, shown in Figure 4B, the motor torque is independent of temperature because 
the interaction potential V is independent of the temperature. For the temperature range 
(AT/T ~ 10%) considered here, ko- and k + are approximately constant. However, both 
&o+ and ko- depend on temperature through a torque-dependent factor: expl—rSi/ksT), 
where Si can be interpreted as the size of the (angular) conformational change of the stator's 
transition state before stepping. A lower bound for Si can be estimated from the maximum 
torque r and the value of r: Si > [/cbT/2t ] ln(l/r), which is roughly 0.2 — 0.5° for r = 0.5 
and r = 175 - 500pN-nm. 

The general properties of the torque(force)-speed relationship and the special case 
of the CW rotation in BFM 

The most important and measurable biophysical characteristics of a motor is its force- 
speed curve (for linear motor) or torque-speed curve (for rotatory motor). What determines 
the overall shape the torque-speed curve? Two important factors are the detailed forms 
of the interaction potential and the stepping rate. Here we study the general relationship 
between the these two key factors and the shape of the torque-speed (r — u>) curve. In the 
case of the V-shaped and the semi-parabolic potentials (V-shaped with a smooth bottom), 
which we have focused on so far, the r — u curve is concave, see for instance Figure 6A. 

Now consider a pure parabolic potential as shown in the inserts of Figure 5A and B and 
the following stepping rates: 



k f (A6i) = kofexpi-rM/keT), 



(12) 



k b (A9i) = k ob exp(TiS9i/k B T), 



(13) 
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which depends exponentially on the torque Tj(= —V'(A9i)) for the forward and backward 
steps with constant k Q b, k f and 69i. Figure 5A shows the torque-speed curves for different 
values of N. The torque-speed curves are slightly convex, almost linear. Another interesting 
finding is that for this choice of the stepping rate dependence on torque, the maximum 
speed at low load increases with N. This dependence of u on N is a consequence of the 
unbounded exponential increase of the forward jumping rate with the negative torque. As 
the number of stator increases, the positions of the stators in the waiting phase can be 
pushed into regimes with larger values of negative torque due to the quadratic nature of the 
potential and henceforth larger jumping rates according to equation ( I12p . Consequently, the 
larger forward jumping rates in the waiting phase lead to higher maximum speeds for larger 
N. 

In order to maintain the independence of uq on N, the jumping rate r(rj) needs to saturate 
to a maximum value as |rj| increases. For example, we can use the following expressions for 
the forward and backward jumping rates: 

k f (A9i) = koflexpifaSet - E )/k B T)/(l + exp^rM - E )/k B T))}, (14) 

h(A6 t ) = ko^expi-irM - E )/k B T)/(l + exp^rM - E )/k B T))} . (15) 

Here, E Q and S6i set the force scale r c = E /59i beyond which the forward jumping rate 
saturates. Choosing r c ps 0.1t the independence of the max speed on the stator number is 
recovered as shown in Figure 5B. This behavior seems to be independent of the specific values 
of r c and 59[. Supplementary Figure S5A shows he case of r c = 0.2r and Supplementary 
Figure S5B the torque-speed curves for different values of 59i. 

Within this framework our model is able to generate a torque-speed curve in quantitative 
agreement with the CW torque-speed curve observed experimentally 18|. Until very recent 
it has been assumed that CCW and CW rotation are symmetric and exhibit the same torques 
and speeds |45j . However, Yuan et al. measured the torque-speed relationship for an E. coli 



strain locked in CW rotation and found a quasi-linear curve. As discussed in 18], the CW 



torque-speed relationship can be explained within our model by having k + close to zero (or 



much smaller than k_) with a V-shaped interaction potential (see Figure S6a of ref. |31|). 
or by having a parabolic stator-rotor potential. Figure 5C shows an excellent agreement 
between the experimental data and our simulation of the torque-speed curve with bounded 
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jumping rates. 

The stepping-rate ratio r not only affects the dependence of the zero-load speed on the 
number of stators, it also changes the characteristics of the torque speed dependence at low 
load. Figure 6A shows the torque-speed curves for two different values of the ratio r: r = 2 
and r = 0.5. For r = 0.5, the torque r generated by the motor decreases as its speed u 
increases. However, for r = 2, the slope of torque-speed curve changes its sign after the 
knee region from negative to positive, giving rise to a non-monotonic dependence of speed 
on torque (or load), (the same behavior is observed when stator springs are added, see 
Supplementary Figure S4A and B). 

The mechanism for this behavior can be understood within our model. As described in 



31] , the dynamics of the motor can be characterized by the two timescales: the average 



moving time < t m > and the average waiting time <t w >. Their dependence on the motor 
speed are shown in Figure 6B for both r = 0.5 and r = 2. 

These two different forms of the torque-speed relationship can be understood by inves- 
tigating the different dependence of (t w ) on the load £l for r = 0.5 and r = 2 as shown 
in the inset of Figure 6B. At a specific value of the load close to the "knee", where 
the torque-speed curves of the two cases (r = 0.5 and r = 2) start to diverge, the depen- 
dence of the average waiting time < t w > on £l also becomes qualitatively different. While 
< t w > (r = 0.5) continues to decrease and eventually level off at roughly 60/is near zero 
load, < t w > (r = 2) starts to increase at to a higher value ~ 110/xs, almost twice the 
value in the r = 0.5 case. This is caused by the fact that for r > 1 it is harder to step 
(switch-hand) when the stator spends more time in the negative torque regime at higher 
load because fc_ < k + . This increase in waiting time for £x < leads to the decrease in the 
motor speed as the load is lowered (below ££) for r = 2 as seen in Figure 6A. In Figure 6C it 
is shown the sign changing in the torque-speed derivative after the "knee" . The torque-speed 
derivative goes from positive to negative when r goes from 1.5 to 0.5, with a critical value 
r c = 1, where the derivative is not defined (oo), that divides the positive from the negative 
derivative case. 

In the model studied in Q' the stepping rate has a complicated dependence on the 
relative angle between the stator and the rotor and in general favors stepping in the positive 
torque regime, which corresponds to r > 1 in our model. Indeed, the torque-speed curve 
shows a change of sign in the torque-speed derivative at around the "knee" , consistent with 
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what we observed for r = 2 in our model, see for example Figure 2B in [27j]. Experimentally, 
such a non-monotonic torque-speed dependence has never been observed. Therefore, any 
realistic model should have r < 1 (or the equivalent), where the FDS mechanism is relevant. 

DISCUSSION 

In this paper, our original mathematical model framework for the rotary flagellar motor 



3l| is extended substantially to incorporate key biological ingredients, such as the stator 
back-stepping, the stator springs, the temperature dependence of the stepping rate, and 
the CW rotation. Our model allows us to learn about the relevant microscopic mechano- 
chemical details of the flagellar motor, such as the power-stroke potential and the stator 
stepping mechanism, from the "macroscopic" measurements, in particular the torque-speed 
dependence. Summary and discussion of specific findings are given below. 

Backward rotation. The introduction of a torque-dependent back-stepping rate repro- 
duces the continuity of the experimental torque-speed curve near stall (high load). We pro- 
posed that the back-stepping probability is dominant in the far-positive side of the V-shaped 
potential, where the forward-step is negligible. Since the landing points of these back-jumps 
are still on the positive side of the potential with the same positive torque r , inclusion of 
back-jumps in our model naturally explains the observed continuous torque-speed relation 
across the stalling point, when the motor is driven backwards. The torque-speed relation- 
ship in the bacterial flagellar motor is linear up to 100Hz in either direction, which rules 



out models 
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that only allow the motor to move backward over a potential bar- 



rier. Back steps were also considered in [27|]. However, since there the back stepping rate is 
only controlled by the pmf it does not prevent the motor from climbing up and sometimes 
crossing the potential barrier. A more detailed discussion of alternative models is given in 
the Supplementary Information. 

Effects of stator springs. Thorough analysis of parameter dependence of speed at near 
zero load shows that adding the stator spring preserves our original model results for all 
ranges of spring constant values as long as r < 1. However, the stator spring constant 
needs to be highly restricted when r > 1. In particular, even a very conservative value for 
the load-free speed variation |1 — u(8)/u(l)\ = 15% would constrain the torsional stator 
spring constant into a very small range. Furthermore, the weak stator spring required for 
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r > 1 gives rise to unrealistic large displacements and fluctuations of stator positions, which 
further limit the validity of the soft stator spring scenario. 

Characteristics of the torque-speed relationship. Performance under external load gener- 
ally characterizes all motor proteins, including both rotary motors and linear motors. In 
our simple model, the torque-speed curve is determined by two factors: the power-stroke 
potential and the dependence of jumping/stepping rates on the mechanical coordinate. For 
example, we show here that changing the ratio of stepping rates r = k+/k- can change 
the torque-speed dependence from monotonic to non-monotonic; and changing the potential 
from linear to quadratic can leads the torque-speed curve to go from concave to convex. 

The newly measured torque-speed curve for the CW rotating BFM [18| can be quanti- 
tatively reproduced by our model with a parabolic stator-rotor interaction potential. As 
shown in our previous work 31(, such a liner torque-speed curve can also be obtained with 



a V-shaped potential by setting the jumping rate k + equal or near to zero. Is it possible 
to distinguish between the two scenarios? In the low to medium load regime where the 
recent experiments have been done, these two mechanisms generate the same quasi-linear 
torque-speed relationship. However, at high load these two scenarios give rise to oppo- 
site torque-speed curvature. Specifically, the parabolic potential shows an increasing of the 
torque much faster than linear whereas the small k+ scenario leads to a constant torque near 
stall. A detailed measurement of the torque-speed near stall for the CW rotation will be 
able to distinguish these two scenarios. The experimental data are fitted better by having a 
bounded jumping rate, which (as shown in Fig. 5B) leads to the independence of the speed 
near zero load on the number of stators for the CW rotating motor, just like for the CCW 
case. This prediction can be checked with resurrection experiments for CW motor at low 
load. 

Finally, we notice that the convex shape of torque-speed curves for the parabolic operating 

B';ential are very similar to the force- velocity curves of linear motor kinesin-1 at high load 
, [l?]], also the collapse at maximum velocity agree with current theoretical model for 



linear kinesin-1 motor [46] and for linear myosin motor systems [471 . These similarities 
suggest that our simple modeling framework may be used to understand/categorize the 
general force-speed dependence for all the mechanical motor molecules. 

Chemical and mechanical gating of the stepping process. Within our model framework, 
temperature affects flagellar motor operation through modifying the stepping rates. Such 
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thermally activated process can be formally treated by the transition state theory [48 ]. 
From the value of the free energy barrier AG obtained from the temperature dependence 
of the stepping rate, stepping seems to be predominately controlled (gated) by the pmf. 
Consequently one would expect that the pmf dependence of the motor speed near zero load 
could have an exponential component. More detailed theoretical models, which consider 
gating mechanisms where proton translocation is modeled explicitly by a barrier crossing 



event, report different behaviors. In 
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271 ] the barrier height depend only on pmf and the 



speed at low load increases almost linearly with the membrane voltage. In 28) both torque 



and speed were predicted to grow sub-linearly with pmf. None of these predictions can be 

inear pmf dependence in a relatively small 



491 ] . We believe direct measurements of the 



ruled out by the indirect experiments showing a 
range of rotation speed between 80Hz and 270Hz 
pmf dependence of the motor behavior at near- zero load are needed to elucidate the exact 
chemical gating mechanism of the stator. Besides pmf our study shows that there must also 
be a mechanical contribution to AG, which accounts for the torque-dependent stepping rate. 
Depending on the value of r, this mechanical gating energy is small (~ knT for r = 0.5) 



in comparison with the pmf, but its effects are significant as shown in [3l[ and here. A 
lower bound for the conformational change (~ 0.2 — 0.5° in angle or ~ 0.8 — 2 angstrom in 
distance) in the transitional state of the stator can be deduced from the mechanical gating 
energy and the maximum torque. 

The pre-factor ko in the stepping rate in equation ( fTTT) . often interpreted as the intrinsic 
molecular fluctuation frequency or the reaction attempting rate, should be proportional to 
\fjjrn, where / is the stiffness of the harmonic energy landscape of the reaction and m is 
the mass of the reacting molecule. For flagellar motor driven by the proton motive force, 
the stepping rates, when keeping temperature constant, should be proportional to 1/y/m. 
Hence, the intrinsic maximum speed of flagellar motor near zero load should be scaled by \/2 
when changing the operating environment from water (H2O) to heavy water (D2O). This 



isotope effect has indeed been confirmed by recent experiments 44]. 
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Figure Legends 



Figure 1. The motor dynamics with stator springs. (A) A a sketch of the model 
for the flagellar rotary motor with the load (blue dot), linked to the rotor (green dot) with 
a spring, and a stator linked to the peptidoglycan cell wall (PG) by a spring. The red line 
represents a possible potential describing the stator-rotor interaction. (B) The stator, rotor 
and load dynamics at high load (N=l, r=2), with ks =200pN-nm-rad~l and D$ =500rad 2 /s. 
The black solid line shows the stepping dynamics of the stator's chemical coordinate, the 
red dashed line represents the stator's overall dynamics: internal coordinate jumps plus 
the physical angular movements. The stretching of the stator spring at steady-state is 
roughly 105. Green and blue dashed lines represent the rotor and load angular coordinates 
respectively. (C) oo r = u(8)/u(l) as a function of k s with D s = 500rad 2 s _1 . Black dots, 
red squares, and green diamonds, correspond to r=0.5, r=l, and r=2 respectively. The 
black dashed line represents the 20% error line: u r — 20%. The light blue semitransparent 
area represents a zero-load speed variation of 10% (as estimated from the experiment). (D) 
Spring displacement A9 S , red squares, and fluctuations Var(0 s ), black dots, normalized by 
5, during the waiting phase as a function of the spring constant at high load for the case 
N=l, r=2. Reasonable values on the order of two 5 or less correspond to values of k s bigger 
than 800pN-nm-rad _1 (black arrow). On the other side u r inside a zero- load speed variation 
of 10% needs k s to be smaller than 20pN-nm-rad" 1 (light blue semitransparent area). 

Figure 2. Phase diagram of the motor behavior in the parameter space. The 

spring constant ks is on the x-axis, and the ratio of stepping rates r on the y-axis. The 
red line with black circles in the main plot shows the parameter choice to achieve 15% 
zero-load speed variation (e r = 15%), model parameters below the red line correspond to 
e r < 15%. The inset shows the standard deviation (red dots) of individual stator at different 
spring constant for the case of N=8, r=0.5, and £ L = 0.002pN-nm-s-rad -2 (low load). The 
allowed space for stator position fluctuation is set as half of the average spacing between 
neighboring stators, and is 22.5° for N = 8 which gives rise to the minimum spring constant 
to be ~ 488pN-nm-rad _1 (green arrow and green triangle in the main plot). The two 
constrains, the zero load speed and the stator position fluctuation, both control the choice 
of parameters together. For r > 1.0, the stator spring constant is highly limited, as shown 
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by the small dark gray region, which diminishes when r > 1.5. Whereas for r < 1.0, the 
allowed stator spring constant, as shown by the green region, has much broader choices with 
essentially no upper limit. 

Figure 3. Dynamics of the motor driven by external forces. (A) The backward 
forcing scheme. The rotor (green dot) moving on the rotor-stator potential (red curve) is 
elastically linked to the load (blue dot), to which a force is applied in direction opposite to 
the force applied by the stator, given by the slope of the potential. (B) Torque-speed curve 
for N=l near stall. Starting from the black dot, the dashed black and red lines represent the 
cases when the applied external torque has respectively the opposite and the same sign as 
the stator generated torque. The inset shows the potential between rotor and stator (dashed 
green line) and the forward and backward rates are represented by the red and black line 
respectively. (C) The forward forcing scheme. The applied force is in the same direction 
as the rotation of the rotor. (D) The normalized motor torque versus speed as the motor 
is driven forward by external forces for N=l, solid black line, and N=8, dashed red line. 
As the external torque drives the motor rotate faster than its natural maximum speed, the 
motor generates net negative torques. The motor torque eventually saturates at a minimum 
with its absolute value roughly the same as the motor's natural maximum positive torque 
for both N=l and N=8. 



Figure 4. The temperature dependence of the torque-speed curves. (A) Speed 
near zero load as a function of temperature. Solid red light shadow line is from numerical 
simulations of our model. The width of the curve corresponds roughly to the error of the 
maximum speed. Black dots are experimental data from (44). Three points (open circles 



with error bar) are data from 
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(B) Torque-speed curves at different temperatures, 
shown in the inset, with stator number from N=l to N=8. V-shaped potential is used and 
the parameters for the stepping rates are: AG = 52k J/mol, ko- = exp(30A). 



Figure 5. The torque-speed curves with parabolic potential. Different number of 
stators, from N=l to N=8, and different stepping rate are used. (A) Unbounded stepping 
rate, shown in the inset, are given by equations ([12]) and (fl3|) with parameter values k r = 0.01 
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and k = 12000s -1 . Black dashed-dotted line: backward jumping rate; red solid line: 
forward jumping rate; green dashed line: parabolic stator-rotor potential. (B) Bounded 
stepping rates, shown in the inset, are given by equations (I14p and ( II 5p with parameter 
values E c = O.ItqSOi, k r = 0.001 and k = 30000s -1 . Parameter values common to (A) 
and (B) are: Ks=10 5 pN-nm-rad _1 , Dg = 500ra<i 2 s _1 , 59i = 0.15. Black dashed-dotted line: 
backward jumping rate; red solid line: forward jumping rate; green dashed line: parabolic 
stator-rotor potential. (C) Normalized torque versus speed for CW rotation. Green lines 
(error bars) are experimental data taken form Gray transparent line simulation of 

our model, with k Q = 32000s -1 and N=8. All other parameters have the same value as in 
Figure 5B. In order to compare the simulations with the experimental published data, the 
normalization of the gray curve has been done at the speed value corresponding to the speed 
value of the first experimental data, i.e. roughly 61 Hz, and not at its maximum value at 
near zero speed. Notice, that the maximum normalized torque value for the experimental 
data is not equal to 1. In fact, the CW data points in [l8| are normalized respect to the 
CCW curve, not shown here. See [18j for details. 



Figure 6. The dependence of torque-speed curve on r. (A) The torque-speed (r — u) 
curves for two different values of r with N — 1. Red squares r=0.5, black dots r=2. (B) 
The average waiting-time < t w > for r=0.5 (black solid line) and r=2 (green dotted-dashed 
line), and the average moving-time < t m > for r=0.5, red dashed solid line, and r=2, blue 
dashed line, over 500 revolutions as a function of the rotational speed for N = 1. The inset 
shows the zoom of the average waiting time < t w > as a function of the the inverse of the 
load, £i, for the two cases: r=2, black dashed line, and r=0.5, red dotted line. (C) The 
torque-speed (r — u) curves for 11 different values of r, from 1.5 to 0.5. The red solid line 
represent the curve corresponding to the critical value r c = 1. 
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Figure 3 
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Figure 4 
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Estimation of the a-helix spring constant 

In bacterial flagella motor, the stator is anchored onto the peptidoglycan 
cell wall through a short segment of a-helix. When the motor rotates and 
the stator interacts with the rotor, the a-helix exerts force that generates 
a torque to anchor the stator around its posi tion. Thi s phys ical process is 
abstracted as a stator spring in some models ( iBai et aL 120091 ). and we have 
systematically studied the effect of this additional elastic element to flagella 
motor operation in the main text of this paper. In this supplemental section, 
we aimed to analyze the magnitude of the stator spring constant, which 
provides important insight to the physical limitation of determining model 
parameters. 

Under most biophysical conditions, a-helix can be regarded as an elastic 
rod. In general, a-helix can bear bending, twisting and stretchi ng deforma- 



tions . And the elasticity of a-helix has been theoretically studied (ISeungho &: Sean 



20051 ). The persistence length l p of an a-helix rod, which determines the stiff- 
ness of a linear object against deformation, has been estimated to be lOOnm. 
In reality, the 3-dimensional deformation of the stator anchoring a-helix could 
be complicated. A reasonable simplification is to assume that the helix only 
bends within the rotor plane. This assumption enables us to estimate the 
lower limit of the stator spring constant, which is also the key point with 
most debates among different theoretical models. 

For an a-helix with length L, we can use the theory of linear elasticity 
to evaluate its deformation energy due to bending. According to the 2- 
dimensional flagella motor geometry shown in Figure [ST| the elastic energy 
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embedded in the a-helix can be written as: 



E a = h B T J^dsl p -A 2 (s) (SI) 

where A is the torsional angle along the bent a- helix and k B T is the ther- 
mal energy. In the limit of small deformation, the torsional angle can be 
approximated using the end displacement of the a- helix d: A « And 



the energy becomes: 



E a (d) « \k B T -L-l p - (^f = \k B T (S2) 



From Figure IST| it is easy to obtain the geometrical relation: d = R ■ 9. So 
the energy contribution from stator spring in the rotor rotational coordinate 



is: 



E a (6) « h B T (S3) 
And the stator spring constant Kg is ; 

ks = dQ2 ~k B T-l p -— (S4) 

A bacterial flagella motor rotor is about R = 20nm in radius, the length of 
the elastic a- helix is about L = 7nm, and under room temperature, lk B T = 
4.2pN-nm, which gives rise to an estimated value for the stator spring: n$ ~ 
SOOpN-nm-rad^ 1 . 
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Stator spring effects at high and low load and 



model parameters study 

In this supplementary section, we explore how the stator spring affects the 
general behavior of a flagella motor. We focus on the three following points: 

• Stator spring effects on flagella motor dynamics at high and low load 
for jumping rate ratio r = 2. 

• Robustness study of load free speed ratio u r as a function of the stator 
spring constant for different values of the stator diffusion constant. 

• Independence of the speed near zero load on the stator number for a 
large interval of values of parameters model, such the stator diffusion 
constant and the the stator spring, when r=0.5. 

Figure IS2l shows typical cases of the stator dynamics at high and low load 
for different number of stators (N) with presence of stator springs. The stator 
spring const ant is chosen to be Ks — 200pN-nm-rad _1 , in line with the value 



used in ref. f Bai et al 



20091 ). We also choose the stator drag coefficient £s 
to make corresp onding di f fusion constant D$ = ksT/^s — SOOracPs" 1 , the 



same as used in( 



Berg's paper (jYuan fc Berg 



Bai et al 



2009J), which is half of the value used in Yuan and 



20081 ). Figure IS21A presents the stator dynamics 
of TV = 1 at high load. The initial angular positions for the load and the 
rotor are set to be the same. At the beginning of the motion, the load is 
stretched for a transient time of ~ 5ms, after which the averaged speed of 
the internal coordinate becomes the same as the averaged stator speed, with 
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a displacement of roughly 105 between the stator position {9\) and the stator 
anchoring point (91), where 5 is size of a FliG unit. The transient time and 
the stator spring stretching both depend on the load and the stator spring 
constant. The results reported in the main text are all computed after the 
system reaches steady state. 

Figure IB2B shows the behavior of single-stator motor at low load (high 
speed) regime. From the simulated trajectory, it is easy to see that the stator 
motor motion exhibits significant fluctuation, although the stator spring only 
get stretched less than 1 ~ 25. In this regime, the load does not show 
noticeable displacement with respect to the rotor. 

Increasing number of stators influences motor dynamics at high and low 
load regimes in different ways. At high load, the qualitative behavior for a 
8-stator motor is the same as for a single-stator motor. Each spring is in 
average stretched ~ 105 in the same direction as the motor rotation. Figure 
IS2l C shows the dynamics of one of the eight stators at high load. All other 
stators behave the same (data not shown). The synchronized property does 
not exist for motor motion under low load. Each stator is stretched differently 
and there exists faster stator, i.e. the stator that jump more frequently, who 
leads the motor of the rotor (as shown in Figure IB21D). Also the directions 
of stators displacements now are different with some stator spring stretched 
in the same direction as the motor rotation and same others in the opposite 
one. In both directions, the stretching between the rotor and the stator's 
chemical coordinate (and physical position) can be as large as many folds of 
5. To keep the stator fluctuation within 28 under the entire motor operating 
spectrum, the stator spring has to be stiffer than SOOpN-nm-rad" 1 , which is 
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4 times l arger than the maximum stator spring constant (200pN-nm-rad l ) 



used in ( Bai et al 



20091 ) to ensure the independence of the load free speed 
on the number of stators. In the main text, Figure 2, we further loosened the 
allowed limit for stator fluctuation, but the obtained stator spring constant 
(~ SOOpN-nm-rad^ 1 ) would still cause significant (> 30%) deviation in the 
load free speeds at different number of stators. 

In Figure [S3| we present the stator spring dependence of the load free 
speed ratio u r = u(8)/u(l) at different stator diffusion constant D$- We 
evaluated the dependence at low diffusivity D$ = 50rad 2 s _1 and high diffu- 
sivity Ds = 5000rad 2 s~ 1 . The dependence is very similar in these two ex- 
treme choices of Ds and is also similar to the intermediate Ds = 500rad 2 s _1 
as shown in Figure IB of the main text result. 

In the case of r = 0.5, all torque-speed curves, with N from 1 to 8, 
collapse near zero load for all choices of stator spring constant k$ crossing 
from 10pN-nm-rad _1 to 10 5 pN-nm-rad _1 , and of stator diffusion constant Ds 
crossing from 50rad 2 s _1 to 5000rad 2 s _1 . This behavior is shown in Figure 
IS4l A and [S4lC. On the contrary, the torque-speed curves collapse only for 
small values of the stator spring constant, e.g. ks = 10pN-nm-rad _1 , when 
setting ratio of jumping rates to be r = 2 (Figure IB4B and IS4D). It is also 
easy to see from Figure IS4T J and ISit D that increasing stator spring constant 
k,s decreases the load free speed of the 8-stator motor, and the magnitude of 
decreasing is proportionally bigger for r = 2 than for r = 0.5. 
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Comment on an alternative models 



In this supplemental section, we will discuss the aspects of the backward 



mot io n and the st ator spring effect for a . 



ences fjXing et al . 



2006 



Bai et al 



2009 



ternative models propo sed in refer 



Van Albada et al 



20091 ) 



When external load exceeds the stall-torque value, stator would be forced 
to move backwards which overpowers the stator-rotqr int eraction. In the 



models presented in refs. (jXing et al 



2006 



Bai et al 



200% (Figure 4B and 



Figure 1 in the referenced papers), sharp peaks in the interaction potential 
were introduced to prevent thermal fluctuations from taking the motor down 
to the back side of the potential field. This setup inevitably results in dis- 
continuity in the torque value for motor backward motion (a small barrier 
for backward rotation), which is inconsistent with the experimental observa- 
tions for E. coli flagellar motor. Moreover, according to this potential setup, 
after the motor moves across the barrier, the required external torque to 
maintain backward motion would sharply decrease, so at near stall torque 
regime, strong fluctuation on the torque values is expected. We are expect- 
i ng the same behavior al so in the model with periodic potential used in ref. 



flVan Albada et al 



20091 ). However, this prediction has not been observed in 
any experiments. 

Another discrepancy between Bai's model and experiments is the stator 
spring effect on motor dynamics. As pointed out in the main text and in 
previous supplemental sections, for the chemical jumping rate ratio r = 2 , 



the st ator spring constant k$ < 200pN-nm-rad 1 (jXing et al 



2006 



Bai et al 



20091 ) gives rise to a 30% deviation in the load free speeds between single- 
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stator and 8-stator motors. Another unphysical consequence from Bai's 
model is the large fluctuation of stator position due to soft spring (Figure 
EH and Figure 2 in the main text). And the fluctuation could be more than 
2 FliG unit size, which could lead to physical contact of neighboring stators 
during rotation. 



Rotation under external torque for N > 1 and 
effects of backward jumping probability 

In this supplemental section, we extended our study of flagellar motor ro- 
tation under external torque from single-stator (N = 1) to multi-stator 
(N > 1). The stator-rotor potential and the jumping rates are defined in 
the main text. 

We first confirmed that the toque-speed curve is still linear with speed 
up to 100Hz in both rotational directions. Figure [S5TA shows an example for 
N = 8. 

We then computed the entire toque-speed relation in normal motor op- 
erative range when the backward jumping probability is introduced. Results 
show that the major characteristic shape of the torque-speed curves remains 
the same, especially the independence of the maximum speed on number of 
stators is preserved (Figure IS5B). 

It is worth me ntioning; that somet imes motors break when large external 



torque is applied ( [Turner et al 



1996). Some motors break catastrophically: 



the cell suddenly spins at higher speed and then stop spinning when the 
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external torque is turned on. However, most motors broke progressively. 
In this case, motor can be partially or completely recovered. A plausible 
explanation for these behaviors is that the catastrophic breaks disconnect 
the drive shaft, whereas progressive breaks disassemble the torque generating 
units, which can be replaced by the cell on time scale of several minutes. 

Force scale and stator conformational change 
in the jumping rate 

The independence of the speed near zero load on the stator number obtained 
with parabolic potential and jumping rates given by equations (14) and (15) 
in the main text does not rely on the specific value of the force scale r c = 
Eq/SOi used in the expression of the jumping rate. Figure S6A shows the 
case with r c ~ 0.2tq. 

A lower bound for S8i, which can be interpreted as the size of the confor- 
mational change of the stator's transition state before stepping, has been es- 
timated to be 0.5° (main text) from the maximum torque To and the jumping 
rate ratio r , which correspond to roughly 0.055. Figure S6B shows torque- 
speed curves of single-stator motor for different values of 5$i, starting from 
the lower value. For 0.055 < 89i < 0.35 all curves fall in the experimental 
observed values of rotational speed and exerting torque. 
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Figure SI: 2-dimensional flagellar motor geometry. The 26 blue dots represent 
the 26 FliG proteins on the rotor of radius R, rotated of an angle 6 at velocity uj. 
The gray rod represent the a-helix of length L, and d its displacement tangential 
to the rotor. PG is the peptidoglycan cell wall. See also Figure 1 in reference 
(jMeacci & TuL 120091 ) 'or a schematic illustration of the rotor-stators spatial ar- 



rangement. 
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Figure S2: The stator, rotor and load dynamics for the case r=2, 
k s =200pN-nm-rad~ 1 , D s =500rad 2 /s. (A) High load with N=l. The red solid 
line shows the chemical internal coordinate of the stator 0\ and the green dashed 
line the dynamics of the corresponding stator: internal chemical coordinate 0\ plus 
the (physical) angular distance of the stator from its stationary anchoring point on 
the cell wall 0±. The stretching of the stator spring is roughly 105. Black and blue 
dashed lines are the rotor and load angular coordinates respectively. (B) Low load 
with N=l. The red solid line shows the chemical internal coordinate of the stator 
0\ and the green dashed line the dynamics of the corresponding stator. Black and 
blue dashed lines are the rotor and load angular coordinates respectively. The rest 
position (0 — 0l)o is set to zero. (C) One of the eight stator (N=8) at high load. 
The red solid line shows the chemical internal coordinate of the stator of one stator 
i, and the black dots the corresponding Of. The stretching is roughly 105. (D) Low 
load with N=8. The red solid line shows the chemical internal coordinate of the 
leading stator and the black dots the dynamics of the corresponding stator. The 
stretching between the two is roughly 105. The solid green line shows the internal 
coordinate of one of the following stators and the blue dashed-dotted line the cor- 
responding stator. Black solid line and black semi-dashed line, indistinguishable 
in the figure, represent the rotor and the load respectively. 
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Figure S3: (A) ui r = w(8)/w(l) as a function of k$ with D$ = hQrad 2 s~ 1 . 
Black dots, red squares, and green diamonds, correspond to r=0.5, r=l, and r=2 
respectively. (B) uj r = lu(8)/lu(1) as a function of k$ with D$ = 5000r ad 2 s^ 1 . 
Black dots, red squares, and green diamonds, correspond to r=0.5, r=l, and r=2 
respectively. 
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Figure S4: (A) Torque-speed curves for four different combinations of ks and 
D,s, when r=0.5 and N=l. Solid black line ks — 10pN-nm- rad" 1 , D s = 5000. 
Dotted red line ks = 10 5 pN-nm-rad~ 1 , Ds = 5000. Dashed green line k$ = 
10pN-nm-rad -1 , Dg = 50. Dashed-dotted blue line ks = 10 5 pN-nm-rad _1 , D$ = 
50. (B) Torque-speed curves for four different combinations of ks and Ds, when 
r=2 and N=l. Solid black line k s = lOpN-nm-rad^ 1 , D s = 5000rad 2 s _1 . Dot- 
ted red line ks = 10 5 pN-nm-rad _1 , Ds = 5000raci 2 s _1 . Dashed green line ks = 
10pN-nm-rad _1 , Ds = 50rac? 2 s _1 . Dashed-dotted blue line ks = 10 5 pN-nm-rad _1 , 
Ds = hQracPs" 1 (C) Torque-speed curves for four different combinations of 
ks and Ds, when r=0.5 and N=8. Dashed black line ks = 10pN-nm-rad _1 , 
Ds = 50rad 2 s _1 . Dotted red line ks = lOpN-nm-rad -1 , Ds = 5000rad 2 s _1 . 
Solid green line ks = 10 5 pN-nm-rad _1 , Ds = SOracPs . Solid blue line 
ks = 10 5 pN-nm-rad _1 , Ds = 5000r ad 2 s^ 1 . (D) Torque-speed curves for four 
different combinations of ks and Ds, when r=2 and N=8. Dashed black line 
ks = lOpN-nm-rad -1 , Dg = 50rad 2 s _1 . Dotted red line ks = 10pN-nm-rad _1 , 
Ds = 5000ra<i 2 s -1 . Solid green line ks = 10 4 pN-nm-rad _1 , Ds = bQrad 2 s~ l . 
Solid blue line k s = 10 5 pN-nm-rad _1 , D s = 5000rad 2 ,s _1 . 



11 



A 




Speed (Hz) 

Figure S5: (A) Torque-speed curve for N=8 near stall. Starting from the black 
dot, the dashed black and red lines represent the cases when the applied external 
torque has respectively the opposite and the same sign as the stator generated 
torque. (B) Torque-speed curves, from N=l to N=8, for the case with backward 
jumping probability. Parameter values as in Figure 3 of the main text. 
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Figure S6: (A) The torque-speed (r — uS) curves for different number of 
stators. N increases, from N=l to N=8, starting from the black curve at 
low torque value to the dark green at high torque value. Parameter values 
are: k = 30000s" 1 , k s = lO^N-nm-rad" 1 , D s = 500rad 2 s"l, k r = 0.01, 
S8i = 0.15, E c = 0.2tq58i, i.e. r c ps 0.2tq, and with parabolic potential. 
(B) The torque-speed (r — u) curves for different values of 59i. N=l, ks = 
lO^N-nm-rad" 1 , D$ = 500rad 2 s"l, k r = 0.01 and with parabolic potential. 
The jumping rates are given by equations (14) and (15) in the main text. 
Solid black line S9i = 0.055, red dashed line S9i = 0.15, green dashed-dotted 
line 59i = 0.25, and blue dotted line 59i = 0.35. 
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